Generalized Lyapunov Exponent and Transmission Statistics in One-dimensional 

Gaussian Correlated Potentials 
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Distribution of the transmission coefficient T of a long system with a correlated Gaussian disorder 
is studied analytically and numerically in terms of the generalized Lyapunov exponent (LE) and 
the cumulants of InT. The effect of the disorder correlations on these quantities is considered in 
weak, moderate and strong disorder for different models of correlation. Scaling relations between 
the cumulants of InT are obtained. The cumulants are treated analytically within the semiclassical 
approximation in strong disorder, and numerically for an arbitrary strength of the disorder. A small 
correlation scale approximation is developed for calculation of the generalized LE in a general corre- 
lated disorder. An essential effect of the disorder correlations on the transmission statistics is found. 
In particular, obtained relations between the cumulants and between them and the generalized 
LE show that, beyond weak disorder, transmission fluctuations and deviation of their distribution 
from the log-normal form (in a long but finite system) are greatly enhanced due to the disorder 
correlations. Parametric dependence of these effects upon the correlation scale is presented. 

PACS numbers: 05.40.-a, 72.15.Rn, 42.25. Dd, 42.25.Bs 



I. INTRODUCTION 

Properties of stochastic systems can be significantly af- 
fected by the degree of correlation of the external noise 
term (see, e.g., Ref. [1] and references therein). As 
an example of such phenomenon, we study statistical 
properties of Anderson localization 2] in a linear one- 
dimensional model 



dx 2 



+ V(x)ip = eip, 



(1) 



with correlated disorder potential V(x). With a suitable 
change of notation, equation ([T]) can describe a station- 
ary problem either for a quantum particle at energy e, 
or for classical scalar electromagnetic or acoustic waves. 
Equation (p} appears also in many other fields of physics. 
For example, with x considered time, Eq. ([1} repre- 
sents a random frequency oscillator, which is a simple 
paradigm of a stochastic dynamical system (see Ref. Q 
for a more detailed discussion of this and other applica- 
tions). In all these instances, an important quantity is 
the Lyapunov exponent (LE) 7 = lmx r _ i . 00 7 (x) , where 

1 ( x ) = lx nl (l^| 2 + W\ 2 ^j an d ^ is a solution of an ini- 
tial value problem. LE is a non-random quantity inde- 
pendent of the specific realization of disorder V (x) [H . 
On the contrary, " local" LE 7 (x) is random and should 
be described statistically by distribution P (7; x). Distri- 
bution P (7; x) can be studied in terms of its cumulants 
(e.g. Ref. [5]) or using the generalized Lyapunov expo- 

nents A<«> = lim MOO ^ In ^ (\^\ 2 + |V>f) ^ where 

(• • • ) denotes average over the disorder realizations (e.g., 
Refs. |],@). 

LEs 7 and 7 (x) are intimately related to the local- 
ization properties of one-dimensional disordered systems. 
The latter was studied extensively and it is rigorously es- 
tablished that almost all eigenstates in one-dimentional 



(ID) disordered systems are exponentially localized un- 
der rather general conditions Q^ll,!! (see also (3|.(Toj| 
and references therein). The inverse localization length 
of the eigenstates, as well as the asymptotic decay rate 
liniL^oo (2L) 1 InT -1 (L) of the transmission coefficient 
T of the system of length L, are equal |4| to the LE 
7. More generally, statistical properties of the quantity 
(2L) 1 InT -1 (L) and of the introduced above local LE 
7 (x = L) become asymptotically equivalent for large L 
(see Sec. Hi]) . Therefore, in what follows, we discuss the 
problem in terms of the transmission coefficient T, while 
the results apply in a broader context of the local LE 7. 

Transmission coefficient is a basic characteristics of the 
wave transport through the non-uniform media. In the 
electronic systems, T is related to the dimensionless con- 
ductance by the Landauer formula [ll[ . Measurement of 
the transmission coefficients is one of the ways to inves- 
tigate different aspects of Anderson localization experi- 
mentally. Some indications of Anderson localization were 
observed with li ght microwaves cold atoms 
and ultrasound [15]. In most of these experiments, the 
correlation scale of the randomness is comparable to the 
scattered wavelength and should be taken into account to 
obtain a full quantitative description. In the experiments 
in Ref. [14J, localization of cold atoms was found by ob- 
serving localized density profiles, rather than transmis- 
sion. Recently, possible experiments on cold atom trans- 
mission through the disordered optical potentials were 
discussed 16]. It is worth noting that in the actively 
developing field of cold atoms, one of the ways to intro- 
duce disorder is by using the laser speckle intensity pat- 
terns, which are highly correlated [17[. On the theoreti- 
cal side, there is an increasing interest in the disordered 
systems with correlated potentials. In particular, in ID 
systems, disorder correlations were found to have strong 
effect both on the localization length [3] . fl8j . (l9j . (2fjj . plf . 
on the generalized LE [22j],[23[ and on the transmission 
statistics [24|], 25] and, in special cases, even lead to the 
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appearance of extended states |26| . 

The main subject of the present work, is the effect of 
disorder correlations on the transmission statistics in the 
asymptotic limit of large L. Transmission coefficient of 
a long disordered system exhibits large sample-to-sample 
fluctuations and is not a self-averaging quantity. Gen- 
eral properties of the transmission distribution have been 
considered in the past using the composition rule for the 
transmission of a one-dimensional chain of statistically 
identical and independent (or only weakly dependent) 
random scatterers [2j|,|28|]. It was shown that a conve- 
nient variable to deal with is InT (L), since it can be rep- 
resented as a sum of independent (or weakly dependent) 
random variables, for which the conditions of the Central 
Limit theorem, or its modification for weakly correlated 
variables, are valid. Then, by additivity of cumulants 
of independent variables, it follows that all cumulants of 
the In T-distribution, P (hiT; L), grow at most linearly 
in the system length L: 



(([lnT- 1 / 2 (L)] n ))=c n L+o(L) 



n = 1,2, 



(2) 



Here ((X n )} denotes nth-order cumulant of the quantity 
X, and the asymptotic cumulant coefficients c„ are con- 
stants, which depend on the microscopic properties of the 
system. The above result applies to systems with both 
uncorrelated and correlated random potentials, provided 
the considered system can be divided into blocks, which 
are much longer than the disorder correlation scale and 
can be treated as nearly independent scatterers. 

Thus, the general form of the large- L limit of the trans- 
mission distribution is well understood, and the remain- 
ing questions are about the dependence of the cumulant 
coefficients c n on the parameters of the problem. Accord- 
ing to Eq. @, transmission statistics can be discussed 
at the following three levels of "resolution". First, there 
exists a self-averaging quantity (2L)~ l InT -1 , which has 
a non-random limit 



lim (2T) 1 \nT~ l 

L—^oo 



ci = 7. 



(3) 



At this level, all the information about the transmission 
fluctuations is lost in the asymptotic limit L — > 00. Next, 

1 /2 

one can define a variable (InT — (InT)) / (4C2-L) , 
where (InT) = 2c\L : for which the Central Limit theo- 
rem holds, and whose cumulants of the order n > 3 vanish 
in the asymptotic limit L — > 00 as L-("" 2 )/ 2 . Therefore, 
this veriable has a limiting Gaussian distribution with 
zero mean and unit variance. Limiting probability dis- 
tributions and their relation to universality and scaling 
have been discussed in the past 28], 29]. In the present 
context of a ID problem, only two model-dependent pa- 
rameters, ci and c 2l are "remembered" in the limiting 
distribution of (InT - (InT)) / (4c 2 L) 1/2 . This fact is re- 
ferred to as " two-parameter scaling" [28] . In special con- 
ditions, namely when the combined scatterers are weak 
and their reflection phases are distributed uniformly, the 



composition rule yields c\ — c 2 up to weak disorder cor- 
rections, which is a manifestation of single parameter 
scaling (SPS) H3, |2j|. 

Finally, for large but finite L, one can study the non- 
Gaussian corrections to the limiting Gaussian distribu- 
tion, which are characterized by the asymptotic cumu- 
lant coefficients c„, n > 3. The latter are responsible for 
the extreme fluctuations of the transmission coefficient. 
These large deviations can also be studied in terms of 
the generalized LEs, which describe asymptotic growth 
rates of the moments of the inverse transmission coeffi- 
cient T _1 and can be expressed as a sum over c n 's 

Under conditions of weak scattering and phase ran- 
domization, in addition to the SPS relation c\ = c 2 , it 
was found in some models of disorder that cumulant ra- 
tios c n /c\, n > 3, vanish up to the weak disorder cor- 
rections [E|.[30|. Vanishing of higher cumulants is also 
indicated by the coincidence of LE 7 and the second or- 
der generalized LE A in the lowest order of the weak dis- 
order expansion, which is valid for a wide variety of ID 
disordered models 0,11, 0,111. Thus, quite generally, 
cumulant coefficients c„ satisfy 



ca/ci = 1, c„/ci = 0, n > 3, 



(4) 



which holds up to the weak disorder corrections. 

Considerable analytical and numerical effort was de- 
voted to verification of SPS and determining lim- 
its of its validity for various models of disorder 
@,H,0,[m,&[il,Hl. Usually, SPS holds for suf- 
ficiently weak disorder, when the localization length is 
much larger than other scales in the problem (see also 
discussion in Ref. [Hj]). 

Deviations from relations ((3]) were studied for uncorre- 
lated disorder 3] , [5] , 35] as well as for some specific mod- 
els of correlated disorder potentials [HJ|, (25|. In particu- 
lar, for the continuous white noise model, it was shown 
[f| that c 2 /ci = 1 for large positive energy e (weak disor- 
der), becomes slightly different from unity near e = and 
vanishes for large negative e, while the "non-Gaussian" 
ratios c n >^/ci have numerical values of about a fraction 
of unity near e = 0, and vanish for large negative or 
positive e (energy scale is determined by strength of dis- 
order). Thus, in the white noise model, unless e is suffi- 
ciently large and negative, the relative width of the Gaus- 
sian bulk of the In T-distribution depends weakly on the 
disorder strength, while finite higher cumulants indicate 
that some non-Gaussian corrections to the distribution 
tails develop in strong disorder. The latter was also ob- 
served in Ref. [3| , where distribution of local Lyapunov 
exponents was studied in terms of the generalized LEs 
for a continuous 5-correlated disorder. 

Presence of disorder correlations can significantly en- 
hance deviation from relations ([J). For instance, one 
finds \c 2 jc\ — 1| ~ 1 near certain energies of the tight- 
binding model with a weak short-range correlated disor- 
der [2J]. Another example of the effect of correlations 
is scaling relation C2/C1 ~ c\R c , obtained in Ref. (25| 
for strong exponentially correlated dichotomous disorder 
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with a correlation scale R c . This scaling, which replaces 
Eq. (U)), is valid when C\R C ^> 1, which means c^Jcx ^> 1, 
so that Gaussian bulk of the distribution P (InT; L) be- 
comes much broader than in weak or white noise disorder. 

These special examples suggest that disorder cor- 
relations can have a significant effect on the trans- 
mission statistics. One can distinguish two limiting 
regimes. In weak disorder, which can be treated per- 
turbatiyelv, LE depends strongly on the correlation scale 
Rc 0, [II, [11, [Iff . However, it usually remains the 
only parameter of the transmission distribution, and rela- 
tions (UJ) stay valid [30] , [33] , [34] . Another special regime, 
which can be called a semiclassical regime of strong dis- 
order, occurs when typical random potential barriers are 
higher than energy e and sufficiently broad, so that tun- 
neling probability already through a single random bar- 
rier becomes small. In this case localization is dominated 
by the under-barrier tunneling, and one can generalize 
the semiclassical approach of Ref. [25] to study the dis- 
tribution of In T (L) in terms of the statistics of disorder 
excursions above the level V (x) = e. In the present work 
we use this approximation for the strong disorder regime 
to show that c n /c\ ~ (c\R c ) n 1 . 

The main focus of the present research is, however, 
on the effect of disorder correlations on the transmission 
distribution at the transition between these two limiting 
regimes. We consider model ([1} with Gaussian disorder 
V (x), whose correlation function is arbitrary, except for 
the assumption that it is parametrized with a finite corre- 
lation scale R c . The transmission statistics is studied in 
terms of the dimensionless ratios cijc\, Cs/ci and A/ci. 
Quantity A/ci characterizes extreme fluctuations of the 
inverse transmission coefficient T^ 1 , since it is related to 
the ratio of the average to the typical values of T -1 . In 
addition, it can be used to show violation of relations Q , 
because A/ci = 1 when these relations hold [3],[22j. 

To calculate A, which can not be treated exactly in a 
general case of correlated disorder, we develop a small- R c 
approximation 36] . A simple regularization is proposed 
to extend this method to the power law correlations, for 
which a standard approximation turns out to be inappli- 
cable. Numerical simulations are used to verify the an- 
alytical approximation for A and to calculate cumulant 
coefficients ci,2,3 for different types of correlations and 
for different values of R c (and e > 0). As expected, the 
obtained results for Cfc/ci, c$/c\ and A/ci show that SPS 
and Eq. Q hold in the weak disorder limit. Beyond this 
regime, we find that disorder correlations strongly en- 
hance deviations from relations @, as compared to the 
case of white noise, and lead to increase of both typical 
and extreme fluctuations of transmission. Dependence 
of these effects on the parameters kR c and c\R c is dis- 
cussed (k = y/e is the wavenumber). In general, effect 
of disorder correlations on the transmission statistics be- 
comes important when c\R c , i.e. the ratio between the 
correlation and the localization lengths, is comparable to 
or greater than unity. 

The outline of the paper is as follows. In Sec. [XT] we 



define the model in question and discuss in some more 
detail the standard and the generalized Lyapunov expo- 
nents. The small-i? c approximation for the generalized 
LE A is developed in Sec. IIIII First, using the analogy 
of Eq. ([1]) with the Langevin equation, equations for the 
second moments of tp (x) are obtained by means of the 
Furutsu-Novikov formula. This leads to an infinite hier- 
archy of the integro-differential equations for the second 
moment of the wave function and its functional deriva- 
tives [H,[36j]. Then, the small- R c closure approximation 
is applied and the result is compared to numerical simu- 
lations. In Sec. IIVI we obtain scaling relations for c n /c\ 
in the semiclassical regime of strong disorder. In Sec. [V] 
analytical and numerical results for A and c n are used to 
discuss properties of the asymptotic distribution of the 
transmission coefficient and deviation from relations Q . 
The results are summarized in Sec. I VII Known results 
on the generalized LE for the (5-correlated disorder and 
on the Born approximation are presented in appendices 
lAl and |Bj respectively. Method of numerical simulation 
is explained in Appendix [Cj Some auxiliary formulae for 
the semiclassical treatment of strong disorder are pre- 
sented in Appendix ID] 



II. BASIC EQUATIONS 

We consider one-dimensional Schrodinger equation 



dx 2 



+ V(x)ip = eijj, 



(5) 



where e is energy and V(x) is zero-mean Gaussian dis- 
order potential. The latter is completely described 
by its two-point correlation function C 2 {x — x') = 
(V(x)V(x')), written as 



c 2 (x - x') = v 2 r 



Rc 



(6) 



where Vq and R c are disorder variance and correlation 
length, and the dimensionless function T(x) decays on 
the scale of unity and is chosen so that 



r(o) = 1, 



+00 



T(x)dx = 1. 



(7) 



The white noise model is obtained by taking the limit 
R c — > 0, while keeping R C V Q 2 constant: 



lim C9 (x — x ) = qS(x — x), 
R c — >o y ' v ' 



(8) 



where we have introduced disorder intensity parameter 
9 ' 



C 2 {x)dx = 2R C V 2 



(9) 



The transmission properties of a one-dimensional system 
of finite length L are described by transfer matrix [37j 



f 



. (10) 
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which relates amplitudes of the incident and the outgo- 
ing waves (a time reversal invariance is implied). Here 
(f> r and 4>t are reflection and transmission phases respec- 
tively. In the following sections, generalized LE A and 
asymptotic cumulant coefficients c„ are studied in terms 
of the amplitude A 2 (x) — \ip (x)\ 2 + \tp' (x)\ 2 , where ip (x) 
is a solution of ([5]) satisfying some generic initial condi- 
tions, e.g. ip (0) = 1, (0) = 0. The obtained results 
are valid also for the transmission coefficient T, since 
in the asymptotic limit L — > oo it becomes statistically 
equivalent to the amplitude A (x). This is because A (x) 
is expressed in terms of the transfer matrix (|f 0|) . which 
asymptotically factorizes into the product of large factor 
T _1 and statistically independent of it matrix of the or- 
der of unity. Such factorization takes place, since for L 
much larger than the localization length, phases <j) r and 
4>t become statistically (almost 38]) independent of T, 
while coefficient T becomes exponentially small in most 
realizations of disorder. Thus, for large L, one can write 



expression 



f 2 n 

A = - V —rC n . 



4 <-~< n\ 

n=l 



(15) 



For large L, value of (T 1 (L)) is dominated by the 
long tail of the T _1 -distribution, which becomes highly 
skewed and heavy-tailed. Therefore, rather than repre- 
senting typical values of T -1 , generalized LE A provides 
useful complementary information about the low-T tail 
of its distribution. To this end, it is convenient to define 
the following quantity 



A 

ci 



1 OO 



2 n r 



4 *> — ' n\ ci 

n—l 



(16) 



In A 2 (L) = InT- 1 (L) + O (L°) , 



(11) 



where the term O (L°j , which absorbs phases 4> r: t and the 
initial conditions for ip (x) , is statistically independent of 
the first one, and its cumulants are of the order of unity. 
Therefore, using additivity of cumulants of independent 
variables, one has (cf. Eq.@) 



whose meaning is twofold. On the one hand, deviation 
of the value of p from unity gives an integral measure 
of deviation from the universal weak disorder relations 
(|4|, since p = 1 when these relations hold (cf. Ref. [3|). 
On the other hand, p can be used to characterizes the 
extreme relative fluctuations of the inverse transmission 
coefficient T _1 in terms of the ratio between the mean 
and the typical values of T~ x : 



= lim 



«ln»T-i» 



2 n L 



lim 



((In- A 2 )) 



2 n L 



(12) 



exp (InT- 1 ) 



_ e (2p-l)2Z^ 



(17) 



A. Equation of motion for moments 



Thus, for large L, inverse transmission coefficient T 1 is IIL GENERALIZED LYAPUNOV EXPONENT A 

statistically equivalent to the amplitude A 2 (x = L), and, 
thus, to the local LE 7 (x — L), defined in the introduc- 
tion. In particular, LE 7 coincides with ci, and in the 
following we will use the latter notation. 

As already mentioned in the introduction, cumulant 
coefficients c„ describe asymptotic form of the distribu- 
tion P(lnT;L). Introducing the dimensionless length 
I = c\L measured in units of the localization length cj -1 , 
cumulants of InT can be rewritten as (cf. Eq. @) 



In™ T- 1 ' 2 (I) 



ci 



1 + 0(1°) 



(13) 



Therefore, it is convenient to discuss transmission dis- 
tributions in terms of the dimensionless ratios c„/ci, 
since meaningful comparison of transmission fluctua- 
tions in different models of disorder should be done 
for the same dimensionless system length I. Then, 
for instance, relative fluctuation of InT is expressed as 
«ln 2 T))/(lnT) 2 = ( C2 / Cl )/L 

Another quantity which characterizes transmission dis- 
tribution is second order generalized LE A defined by [6] 

A = lim -L In (A 2 (L)) = lim -L In (T _1 (L)) , (14) 

Since \a.(A q (x)) = In (e qlnA ( x ^ is a cumulant generat- 
ing function of the distribution P (In A; x) , the general- 
ized LE is related to all cumulant coefficients c n by the 



The linearity of the Schrodinger equation ([5]) allows 
one to obtain closed-form equations for the 2n-order 
products [tp(x)] k [tp' (x)] , where k + l = 2n, k,l = 
0,1,2,.... To this end, we rewrite Eq. in the form 
of the Langevin equation. The x coordinate is consid- 
ered as a formal time on the half axis x = t G [0, 00) 
and the dynamical variables u{t) = ip (x) , v(t) = ip' (x) 
are introduced. In these variables the Langevin equation 
reads 



d t u = v, d t v = [V(t) - e]u , 



(18) 



where V(t) is now the correlated noise. We need to con- 
sider an initial value problem for the second order mo- 
ments ^|it| 2 ^}, v / 2 _1 (u*v + uv*), ^|w| 2 ^), whose asymp- 
totic exponential growth rate gives the generalized LE 
A, Eq. (fl"4"]) . Introducing vector 



\u v + uv 



\ 2 \T 



(19) 



and using Eq. (fTS)) . one obtains 

d t Y = (C + V (t) V) Y, Y (0) = Y , (20) 



5 



where 





' 10" 




"0 





0" 


C = y/2 


-e 1 


; 23 = 


1 










-e 







1 






(21) 



and Y is an initial condition. To obtain equation for 
(Y(t)), Eq. ([2U)l is averaged over the disorder, which 
introduces correlator (V (t) Y (t)): 

9 t <Y)=C<Y)+23<y( i )Y>, (Y(0))=Y . (22) 

Applying the Furutsu-Novikov formula [J, [HI yields 

d t (Y)=C(Y)+ J*C 2 (t- T )v/ 5Y {U 



SV (r) 



dr, (23) 



where gyfy is a functional derivative of Y (t) with re- 
spect to V (t) and C2 (£) is the disorder correlation func- 
tion |6]). Since equation (|23|) contains a new quantity 

g v ( * ) ) 1 ^ is not closed with respect to (Y). An impor- 
tant exception is the case of the uncorrelated disorder, 
which is considered in Appendix [AJ and will be used in 
the forthcoming analysis. In general case, to proceed fur- 
ther, an additional equation for gyfy is required, which 
is obtained by differentiating Eq. (f20| functionally with 
respect to V (t) , and using the differentiation property 

Wtfm Y (*) = ftW$) & Thi s yields 
d 8Y (t) 



dt SV (r) 



= VY (t) S (t - r) 



9 (t - r) [C + V (t) V] 



5Y (t) SY (t < t) 



SV (r) ' SV (r) 



= 0, 
(24) 



where O (t — r) is the Heaviside step function. For t > t, 
Eq. can be rewritten as 

95Y ^-[ C + 7 W 23]|XM ) ^) = ^ Y(r) . 



<9i 5V (t) 



<5F (r) <5V (r) 

(25) 

was used in deriva- 



The causality property <5y(r>t) 
tion of Eq. (j2"4l . It follows from the fact that Y (t) obeys 



first order differential equation (|20[) with an initial con- 
dition at t — and, thus, is independent of V (r) at the 
time r later than i. Then, integration of the differential 
equation (|24[) over £ from some t\ < t up to r + yields 
the initial condition in Eq. (|25|) . 



Note that functional derivative ^77^ in Eq. (|25p and 

vector Y (i) in Eq. (|20|) obey the same differential equa- 
tion, but with different initial condition. Due to this 



similarity, moments 



/ «Y(t) 

\8V(t 



and (Y (t)} grow with the 



same asymptotic rate for large (t — r) and t respectively. 
The difference between the two quantities is that the ini- 
tial condition for gyjy^ depends on Y (r) and, therefore, 
it is correlated to V (t). 



Averaging Eq. (1251) over the disorder and using the 
Furutsu-Novikov formula, one obtains the following equa- 
tion for the moment ( 

\ SV(t) 



d_ 

dt 



SV (r) 



_ c /5Y{t) 
\5V(r) 



f C 2 (t- t') V 
Jo 



<5 2 Y it) 



SV (t') SV (r) 



with the initial condition 



SY (t = r H 



= 23<Y(r)) 



dr' 
(26) 

(27) 



SV (r) 

Thus, Eq. (I2"f?l) is again not closed with respect to 
because it contains an additional quantity 



SV(t) 

sv\t^)sv(t) )• Then, the above procedure can be iter- 
ated repeatedly to obtain an infinite hierarchy of cou- 
pled differential equations involving moments of hi ghe r 
functional derivatives of Y with respect to V (t) [lj|.[36l|. 



B. Small-_R C decoupling approximation 

To make progress with the infinite hierarchy of coupled 
differential equation, the first two of which are equations 
(f2"3"|) and (|2"B1) , one usually resorts to some closure approx- 
imation [l|,[36j. In this work we apply a small- R c approx- 
imation to close the second equation (|26p with respect to 

( sv(t) ) ' This allows to express ^ gy(-r) ) i R t erms °f the 
initial condition T> (Y (r)) . Then, substituting it into Eq. 
(f2"3")l . one obtains a closed integro-differential equation for 
(Y). 

A formal solution for gyfy is obtained from Eq. (|25|) 
in the form 



5Y(t) 



SV (t 

where the evolution operator 



K{t,r)VY (r), 



(28) 



/C (i, t) = T exp J dr(C + V(T)V) (29) 



is the T-ordered exponent. Disorder average of Eq. 
introduces a correlator (JC (t, r) 23Y (r)). Application of 
the decoupling, or mean field, approximation yields 



5Y(t) 
SV (r) 



(/C(i-r))23(Y(r)) 



(30) 



where we can write (JC (t, r)) = (JC (t — r)) due to the sta- 
tionarity of V (t). Terms, neglected in Eq. (|3"0)) are given 
formally by the generalization of the Furutsu-Novikov 
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formula for the correlation between two functionals of 
the random Gaussian process [l[: 



(F{V}G{V}) = (F{V}) (G{V}) + 
^ 1 f f I S n F {V} 



n=l 



SV (t n )---SV (h) 



S n G{V} 



Y[C 2 (U - dtidsi. (31) 



SV (s n ) ■ ■ ■ 5V { Sl ) / ±i 



In general case, when functionals ^{1/} and G{1^} de- 
pend on the noise at simultaneous times, the decoupling 
approximation does not rely on small correlation radius, 
but depends on other parameters of the problem, e.g., 
small noise intensity The present case is special, be- 
cause Y (r) and K (t, t) depend on V (t') at separate time 
intervals, < t' < t and r < t' < t respectively. There- 
fore, the decoupling approximation for (fC (t,r) VY (r)) 
is also justified for sufficiently small R c . In particular, 
equation (|30[) becomes exact in the white noise limit (cf. 
Eq. (|37|) below). A small- R c perturbative expansion of 
the n = 1 term in Eq. (|31[) shows that it can be neglected 
if both kR c and AR C are small compared to unity. 

Substitution of the decoupling approximation (|30j) into 
Eq. 423|) yields 

^{Y(i))=C<Y(t)) + 

+V [ C 2 {t- t) {K (t - r)> V (Y (t)) dr, (32) 
Jo 



Note that inserting (Y (t)) = (K. (*)) Y into Eq. 
one arrives at closed nonlinear equation for (K, (t)). Al- 
ternatively, assuming in Eq. (|32[) some explicit approxi- 
mation for (K,(t — t)), which will be defined later, one 
obtains closed linear equation for (Y(t)). It can be 
solved in a standard way by the Laplace transform, 
while its large time it >• R c ) eigenvalues can be de- 
termined by substituting the asymptotic solution in the 
form (Y (t)) = Y^e 4 ^, where Y m is a constant vector 
and A is an eigenvalue to be found. For t 3> R c , owing to 
the decay of the correlation function C2 (t — r) , the lower 
limit of integration in Eq. (|32|) can be replaced with —00. 
Then, changing the integration variable as r — > (t — r) , 
one obtains a stationary eigenvalue problem, from which 
the generalized LE A is found as the largest real root of 
the characteristic equation 



det C + Vi^j C 2 (t) (K (t)> e~ 4AT dTj V - 4A =0. 

(33) 

Equation (|33[) makes sense only when g = 2 J C'2(x)dx 
is finite, which is satisfied by assumption ([7|). This condi- 
tion is important, since the correct asymptotic growth of 
(AC (t)) is proportional to e 4Ar (see note after Eq. ([231) ). 
Therefore, for large r, it should be exactly canceled by 
the factor e~ 4Ar in Eq. (1331) . when A is equal to A. 



White noise approximation for / |^py 



For (K. (t)) in Eq. (|33l> we apply a standard white 
noise approximation [36!] by replacing the true correlation 

function C2 (i — t') in equation f|26[) for ( ) with its 



white noise limit gS (t — t'), where g = 2R C V Q 2 , cf. Eq. 
([!]). The integral in Eq. (f26f becomes 



C 2 (t - r')X> 



5 2 Y (t) 



5V (r') (r) 
(5Y (i) x '* 



SV (r) 



C 2 {t-T')dT', (34) 



A- 2 Y(t 4 



P I y I * l was used 



where the initial condition jymsyjr) 
(cf. Eq. (|25jl). Approximation (|34[) is justified when 
a'V(t')<5V(t) / changes slowly on the scale of R c . Since 
relevant characteristic scales of the problem are wave 
number k — \f\e\ and the generalized LE A, the con- 
ditions to be fulfilled arc 



kR c < 1 and AR C < 1, 



(35) 



i.e. the same conditions as for the decoupling approx- 
imation (I30p . Inserting the white noise approximation 
(|34|) into Eq. ([26)) yields the following closed equation 

for ' 



\<5v(t 



9 /5Y(t) 



dt \ SV (r) 
'SY (r+) 

(jy (t) 



2><Y(t)>, 



*Y(t) 
51/ (t) 



(36) 



which coincides with Eq. (|A1[) . obtained for (Y (i)) in 
the white noise model (Appendix |A"| . The solution of 
Eq. (J36J) is 



*Y(i) 
£F (r) 



= (£™„ (t - r)) V (Y (r)) , (37) 



where {K wn (t - r)) = e ^(^ff)(*-^) an d matrix 



Al (e,^) = C + §X> 2 is given in Eq. (|A2]l . Thus, de 



coupling (1301) is obtained automatically within the white 
noise approximation, whereas the mean propagator in 
(|30| is approximated by the exponent 



(AC (*)) « (AC W » (t)) = e^)*. 



(38) 



For weak disorder, using approximation A4 (e, g) « C in 
the above expression for (jC wn (t)}, and substituting it 
into Eq. (f3"3"|) , one recovers the Born approximation (|B3|) , 
which is valid for an arbitrary R c . It follows that validity 
condition kR c <C 1 is relaxed in the limit of weak disorder 
(while the second one, AR C <C 1, is usually fulfilled). 
Substitution of Eq. ()38j) into the eigenvalue equation 
yields an integral J °° C 2 (r) e M ^ T e~ AAT dT, which 
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is calculated by diagonalizing matrix M. (e, g). The latter 
is not a normal matrix {AiJvO ^ M^Ai), and can be 
written as 139] 



QlM (e,g) Q R =M{e,g)=A 



Ai 
A 2 
A 3 



. (39) 



where eigenvalues A, = 
[5l) , while matrices Ql 



Ai (e, g) are given in Eq. 



i 2 , ix 3 



and Qi 



(u^,uf,u^) = Q^ 1 are composed of the left and the 
right eigenvectors of Ai (e,g), given in Eq. (|A6I) . Insert- 
ing the white noise approximation (|38p into Eq. (1331) and 
using parametrization of the correlation function, Eqs. 
© and one obtains 



det 



C + ^DQrHQlV - 4A 



= 0, 



where the diagonal matrix 



H = h[R c AA-Ai (e,g) 



(40) 



(41) 



is expressed in terms of the Laplace transform of the 
dimensionless correlation function 



h(z) 



/ r(*)e- 

Jo 



s ds, z € 



(42) 



Note that, if the dimensionless correlation function V (s) 
decays slower than exponentially, then h (z) is defined 
only for Rez > 0. As a consequence, the small- R c ap- 
proximation will turn out to be inapplicable in the case 
of slowly-decaying (e.g. by power-law) correlation func- 
tions. This case should be treated with some care and 
will be discussed later in this section. 

From now on, we restrict ourselves to the case e > 
— §<? 2 / 3 . In this regime, eigenvalues of Ai(e,g) satisfy 
A3 = AJ (cf. (|A5p ). and matrix H. can be written as 



u 



hi 
h r + ihi 
h r — ihi 



(43) 



where 



hi=h{ 4R, 



A-Ai 



hi =lmh{ AR C 



Re h AR, 



A-Ac 



A-A; 



(44) 



and h (z) is defined in Eq. (|42[) . Then, the matrix prod- 
uct in Eq. (14TJ1) becomes 



12A{ + e 





fl2 /22 
/l3 fl2 



(45) 



Correlation 


T(x) 


Hz) 


Exponential 
Gaussian 
"Speckle" 


e~ x 

£ -nx 2 /4 

sine 2 § 


e * 2 /T erfc -4= 
2- arctan £ - In 2-^-. Re 2 > 



TABLE I: Dimesionless correlation function L (x) and its 
Laplace transform h(z), Eq. 1)42 p. for the models of disor- 
der implemented in simulations. 



where 



/ 12 = 2 3 / 2 A! (hi - hr) + V2 bAl J e hj, 
3Ai , 



/22 = hi- h r 



/13 = 8A 2 (/ii + 2/i r ) + 2e/i r 



2Aieft-i 

7W+ 



(46) 



Finally, inserting matrix (|45[) into the eigenvalue equa- 
tion (1301) 



det 



C + ^VQrHQlV — 4A 



= 0, 



(47) 



one calculates the generalized LE A as the largest real- 
valued root of the determinant. 

The roots of the corresponding non-linear equation can 
only be found numerically. Simple analytical expressions, 
however, can be obtained in some limiting cases. In weak 
disorder, A 2 <C e, one recovers the Born approximation 



(|B3j) . as already discussed after Eq. (|38|) . In sufficiently 
strong disorder, such that A 2 3> e, one can discard e 
in Eqs. ([46]) . which simplifies the expressions. The re- 
sulting eigenvalue equation can be solved analytically in 
the limits Aii? c <C 1 and Aii? c 3> 1. The first condition, 
Ai_R c <C 1, corresponds to the white noise limit and yields 
A = Ai, as expected. In the second limit, Aii? c ^> 1, 
which is beyond validity of the small- R c approximation 
and is considered formally, one obtains A ss Ai (more 
specifically 1 < A/Ai < 1.16). Between these two lim- 
its of small and large AxR c , one still expects A a A l5 
which was also found solving the eigenvalue equation nu- 
merically for the correlation functions listed in Table HJ 
Thus, in strong disorder, A 2 3> e, we obtain A Ai irre- 
spectively of the form of the disorder correlation function. 
This result was also obtained in semiclassical approxima- 
tion and in numerical simulations discussed below. Ac- 
cording to Eqs. (|A4|) and ([9]), the above strong disorder 
condition can be reexpressed as 



|e| 3/2 « V 2 R C . 



D. A regularized white noise approximation 

In what follows, we refer to the approximation method 
described in the previous subsection, Eqs. (|45|) - (|47)) . as 
the standard white noise approximation (SWNA). As al- 
ready noted after Eq. (|4"0|) , this method is not applica- 
ble to slowly decaying correlations. We now discuss this 
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point in some more detail and propose a simple regular- 
ization to extend the applicability of the approximation. 

The white noise propagator (K wn (t)), involved in ap- 
proximation ([37)) for ^ gy(*) V has the largest eigenvalue 
e 4Ai(e,g)^ ^gj-g Ai (e, 5) is the generalized LE for the 
uncorrelated disorder (see Appendix [A"|) . As follows from 
the Born approximation (|B3[) . it is usually larger than 
the generalized LE A for the correlated potential with 
the same intensity g, at least for sufficiently weak dis- 
order. Therefore, if correlation function C 2 (r) decays 
slower than e _4 ( Al ( e ' 9 ) _A ) r , then the integral in the eigen- 
value equation (|33|) . which can be estimated as 

00 00 

J C 2 (r) {K wn (<)) er iKT dT ~ J ' C 2 (r) e< Al ^)-A)r d ^ 


(48) 

diverges for A = A. In the final expressions, Eqs. 
(l4T))) - (|46p . this potential divergence resides in function 

hi = J °° T (r) e 4 ( Al(c ' s) " A ) r dr, defined by Eqs. 02J and 
flU} . Thus, for slowly, e.g. sub-exponentially, decay- 
ing correlations, A can not be a solution of the eigen- 
value equation and the standard white noise approxima- 
tion is a priori inapplicable. Divergence of the integral 
in Eq. (|48l) as an artifact of the white noise approx- 
imation (1371) . which yields a wrong large-time asymp- 
totics (f$£y) ~ e 4Ai(6, 9 )(t-T) < In deed, as noted after 
Eq. ([25|) . the true asymptotics of the moment of the 
functional derivative is ( sv{t) ) ~ e 4A ( t_r ), i-e. the same 
as for (Y (t)). Therefore, if (JC wn (t)) in (@S) had a "cor- 
rect" asymptotics, then the integral would converge at 
A = A as long as J °° C 2 (r) dr = | < 00. Thus, some 
kind of a regularization is required for t 3> i? c . One 
possibility is to introduce an upper cutoff of the order 
of R c into the inte gra l for hi. Such cutoff was effec- 
tively used in Ref. [23[ for another scheme of a closure 
of Eq. (|23|) . A more simple regularization is just to set 
hi = J Q r (t) dr = 1, which amounts to a self-consistent 
modification of the white noise propagator (IC wn (t)) in 
the approximation (|38[) by replacing its largest eigenvalue 
e 4Ai(e,a)t w i^i! e 4A *. We use the second approach and re- 
fer to it as a "regularized white noise approximation" 
(RWNA). Thus, the regularized white noise approxima- 
tion is given by Eqs. (|45|) - (|471) . with hi replaced by unity. 
Note, that this method can be regarded as applying the 
decoupling approximation (|3T)|) with the corresponding 
effective propagator, different from that in Eq. (|38[) . 

E. Numerical test of small-i? c approximation 

The standard and the regularized white noise approxi- 
mations are compared to the direct numerical simulation 
for the three types of disorder correlations given in Table 
HJ Analytical values of A are obtained by finding nu- 
merically the largest real- valued root of the determinant 




0.01 0.1 1 10 

p 

FIG. 1: (color online) The standard (SWNA) and the reg- 
ularized (RWNA) white noise approximations for the gener- 
alized LE A are compared to numerical simulation for (a) 
Gaussian and (b-d) exponential correlation functions (Table 
[I]). Legends in (a) apply to all panels. GLE, normalized to 
the wavevector k, is plotted against the dimensionless disorder 
intensity j3, Eq. (|49p . Error bars indicate uncertainty of the 
Monte-Carlo calculation of A. For comparison, we also plot 
2 nd order (i.e., one order beyond the Born approx.) pertur- 
bative calculation of A due to Ref. [22|, and the exact GLE 
A wn for ^-correlated disorder given by Ai in Eq. (|A5[l . In 
addition, LE ci, obtained by numerical simulation, is plotted 
with triangles. As expected for Gaussian disorder ;20], A and 
ci coincide in weak disorder as long as 2 nd order perturbation 
theory is valid. In strong disorder, ft ~ 1, GLE for correlated 
disorder coincides with A mn (/3). 



(|40l) . The Monte- Carlo simulations were done for the 
tight-binding model (jCll) with energy E = — 2cos/c near 
the band edge, which translates to the continuous model 
([5]) with a positive energy e = k 2 . Further details of the 
numerical simulations are given in Appendix [Cj 

First, we examine models of disorder with the quickly 
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decaying exponential and Gaussian correlation functions 
(Table IJ). The results are compared to the Monte-Carlo 
simulation in Fig. [TJ The error bars indicate uncertainty 
of the Monte-Carlo calculation of the generalized LE, 
which is explained in Appendix [Cj the generalized LE 
A is plotted for different values of kR c as a function of 
the dimensionless noise intensity 



gk_ V 2 kR c 



Ae 2 



(49) 



where g is defined in Eq. ©. Note that /3 = A/k in the 
Born approximation (|B3I) for the white noise disorder. 

The standard and the regularized white noise approx- 
imations yield close results for kR c < 1, while Difference 
between them increases with R c and, e.g., for the ex- 
ponential correlation with kR c = 3.4, it becomes of the 
order of ten (Fig. [TJd)). Comparison with the simu- 
lation results shows that both approximations are good 
for kR c < 1 (Fig. QJa-c)), but deteriorate with further 
increase of kR c , as expected. For example, for the ex- 
ponential correlation with kR c = 3.4, the resul of the 
standard approximation and the simulated values of A 
differ already by a factor of ten at moderately strong dis- 
order, f3 ~ 0.1 (Fig. Hid)). The regularized white noise 
approximation seems to be somewhat better, though the 
improvement is inconclusive. 

The main advantage of the regularized white noise ap- 
proximation is its applicability to the long-range corre- 
lations, as demonstrated in Fig. [5] As an example, we 
use "speckle" correlation function T (x) = sine 2 | , where 
sine x = sm * x , which describes correlation of the intcn- 
sity in some quasi-one-dimensional laser speckle patterns 
140 1. The latter are used to create disordered poten- 

Note, 



tials in experiments with cold atoms [17] . Note, how- 
ever, that our Gaussian model does not describe a true 
speckle intensity, which exhibits highly non-Gaussian 
fluctuations. This correlation function has a remark- 
able property that the corresponding power spectrum, 
given by the Fourier transform of the correlation function 
2Vq dxT (x / R c ) cos qx , is proportional to the "tent" 
function (1 — \qR c /n\) 6 (1 — \qR c /n\) and vanishes for 
\qR c \ > 7r. As a results, the Born approximation for 
both LE ci and the generalized LE A, Eq. (|B3I) . van- 
ishes for 2kR c > 7r. In view of this fact, the regularized 
white noise approximation is tested below and above this 
threshold for 2kR c equal to 0.97T and 1.03-7T respectively. 
In both cases we find agreement with the simulation re- 
sults within the factor of the order of unity (Fig. 

Let us summarize the above comparison of the ana- 
lytical and the numerical results. The suggested sim- 
ple regularization allows extension of the standard white 
noise closure approximation to the cases of slowly de- 
caying correlations. In the regime of moderate disorder, 
0.1 ^ /3 < 1, white noise approximation is applicable 
for kR c < 1 (we consider e = k 2 > 0), in accord with 
the formal requirement kR c <C 1 (cf. Eq. (f55|) V In the 
limit of weak disorder, j3 <C 1, our approximation coin- 
cides with the perturbation theory (IB3|) . i.e. becomes 




o A, simulation - 

A, RWNA 

■ ' 2 nd order PT 

A ci, simulation i 



0.01 



0.1 
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FIG. 2: (color online) The regularized white noise approxima- 
tion (RWNA) and the numerical simulation of A are compared 
for the power law "speckle" correlation (See Table [TJ> with (a) 
kR c = 0.9tt/2 and (b) kR c = 1.03tt/2 (legends in (b) apply 
to both panels). Error bars indicate uncertainty of the nu- 
merical calculation of A. Like in Fig. [TJ A wn (J3 ) an d the 2 nd 
order perturbative calculation of A due to Ref. [23 ] are given 
for reference. Triangles represent numerical result for LE ci. 
As in Fig. [TJ A coincides with ci in weak disorder, and with 
A„„ (P) in strong disorder. 



exact up to weak disorder corrections. Thus, for /3 <C 1, 
the validity condition kR c < 1 is relaxed (cf. discus- 
sion after Eq. ([35])). Finally, for strong disorder, j3 > 1, 
numerical data coincide with and, thus, confirm the an- 
alytical result A(/3) rs Ai (/3) = A wn ((3), obtained for 
|e| 3/2 < V$Rc, i.e. /3 » 1, at the end of Sec. ImCl This 
is inspite of the fact that for /3 > 1 in our simulations we 
have AR C > 1 , which violates the white noise approxima- 
tion validity condition AR C -C 1, Eq. (|35[) . In the next 
section, however, the same relation A (/?) sw Ai (f3) is ob- 
tained in the regime AR C 3> 1 specifically for Gaussian 
disorder. 



IV. SEMICLASSICAL REGIME OF STRONG 
DISORDER 

In this section we analytically calculate cumulant coef- 
ficients c„, defined by Eq. (|12|) . in the regime when con- 
ditions e < Vq and R c y/Vo — e 3> 1 are fulfilled, where 
Vq is the noise variance introduced in Eq. ([6]). An op- 
posite regime, R c \/Vo — e <C 1, corresponds to the white 
noise limit, for which the cumulants were calculated in 
Ref. Q. We assume that Vq and R c are typical height 
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and width of the random barriers, therefore probability 
of a tunneling through a single typical barrier is expo- 
nentially small under the above conditions. In this case, 
as argued in Rcf. 25], transport is dominated by the 
under-barrier tunneling and interference effects are neg- 
ligible. Since R c is a spatial scale of the disorder V (x), 
the second condition, Rc^/Vq — e 3> 1, justifies the semi- 
classical approximation for the under-barrier tunneling. 
Thus, an exponential growth rate of the solution ampli- 
tude A (x) = \J i\) 2 + ip' 2 can be estimated as a magnitude 
of the imaginary part of the semiclassical action: 



\n\A(L)\ = [ K(x)dx. 
Jo 



(50) 



Here k (x) = e (V) yV (x) — e is an imaginary part of 
the complex momentum, where £ (V) = (V — e) is 
a short notation for the unit step function. Function 
C (V) selects regions where the energy is below the disor- 
der barriers and solution ip (x) grows exponentially, while 
the remaining regions, where ij) (x) has an oscillatory be- 
havior, are discarded. 

Semiclassical approximation (|50l) directly relates dis- 
tribution of In | A (L)\ to the statistical properties of the 
above-threshold excursions of the random process V (x) . 
Disorder average of Eq. ([501) yields LE or the first cumu- 
lant coefficient 



quite involved. It is, however, possible to obtain a simple 
estimate of the multiple integral in (|53|) . As an exam- 
ple, consider the two-point correlator of the regular step 
function (V (x)), which can be calculated analytically 
pl| , and is given by 



((0(y(x))0(y(x')))) 
((© 2 0O» 



— arcsm 

7T 



where the dimensionless correlation function T (x) was 
defined in Eq. (J6)). As expected, the two-point corre- 
lator of 0(V(x)) decays in the same manner as F(x). 
Due to the normalization conditions (0, the integral of 
the expression in Eq. (|54p is of the order of R c . Simi- 
larly, any n-order cumulant in (|53[) has maximum value 
{{n n (0))) and decays with a distance from the origin on 
the scale R c . Therefore, rescaling the integration vari- 
ables in (|53[) by R c , the remaining dimensionless inte- 
gral can be grossly estimated as a volume of the (n — 1)- 
dimensional unit sphere S n -i- This gives the following 
estimate 



S^R?- 1 ««")) 



(55) 



(k) = (n) e <0 e (V)) , 



(51) 



where (n) e = J e dVy/V^eP e (V) is a conditional 
average for V > e over the distribution P e (V) = 
C (V) P (V) I / e °° dVP (V). Here P (V) is distribution 
of the disorder values at a single position. 

Using semiclassical approximation (1501) . higher cumu- 
lant coefficients c„ can be expressed in terms of the cor- 
responding joint cumulants of n (x): 



lim — 



(52) 



i=l 



Since V (x) is a stationary process, the joint cumulants 
((■■•)) in Eq. (j5"2"j) depend only on the coordinate dif- 
ference. The value of the cumulants is maximal when 
all the coordinates coincide, and decays on the scale R c 
with a distance between the points. This is because R c 
is the only spatial scale in our model ©• Thus, shifting 
all the coordinates in the cumulant by x n and assuming 
L » JJ C , the integral in (|51l) can be approximated as 



°° °?n-l / / n-1 

/ •■■ / Hd Xi u K (o)i[ K (x. 

»=1 \ \ i=l 



(53) 



where we have neglected the boundary effects at the cor- 
ners Xi = X2 = ... = x n — 0,L and extended the in- 
tegration to ±oo. For the Gaussian process V(x), an 
analytical calculation of the joint cumulants of k (x) is 



where functions B n (e/Vo) are of the order of unity (ac- 
cording to Eq. dSTJ, B x ee 1), and ((k")> = ((n n (x))), 
which is independent of the position x. Coefficients 
B n >i (c/Vb) compensate for the approximation of the 
(n — l)-dimensional integral by the volume of sphere and 
depend weakly on the specific form of the correlation 
function of V (x) and on the ratio e/Vb- Numerical values 
of B n (e/Vb) can be found in computer simulation by cal- 
culating statistics of the quantity on the right hand side 
of Eq. ([5H|) . and fitting it to Eq. (l55"j) . As an example, 
we obtain B 2 (0) « 0.9, B 2 (|) w 0.8 and B 3 (0) « 1.6, 
B3 (i) w 1.2 for all three types of correlations given in 
Table H 

Combining Eqs. ([ST]) and (|5"5")l . one obtains the follow- 
ing relations 



(56) 



where 



= ((«")) 

- <K.)" 



are dimensionless functions, 



whose specific form depends only on the disorder dis- 
tribution P(V). This result should be contrasted with 
the weak disorder relations ((4]). 

Using the fact that e (V) is either or 1, it is conve- 
nient to express cumulants of k (x) in terms of ((0™ (V))) 
and ((n n )) e , where ((•••))<: denotes a "conditional cumu- 
lant" , calculated with distribution P e (V) [see Appendix 
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D 



For example, for n = 2, 3 one obtains 

«GJ(V)» , 1 «. 2 )) £ 



h TF 



h 77 



(6 £ (V)) (e« (y)) 

(e e (^)) 2 (k)1 
<e e (F)> 3 (k)1 



(Moo)) 

(e e (v)) 3 



(57) 



This expansion is helpful, since it shows separate contri- 
butions of the fluctuation of the indicator function 9 e (V) 
and of the barrier height fluctuation above the level e. 
For the Gaussian distribution, P (V) = 

— 1/ 22/2 / \ 

(2itVq) e~ v ' 2V ° , explicit expression for ci, f% ( ^ J 

and / 3 are given in Eqs. (|D3j) and (|D4|) . 

Finally, the semiclassical approximation ([50)1 can be 
used to calculate the generalized LE A in strong disorder 
limit 



A sc = lim 

L— >ac 



4L 









^exp 


2 / k (x) dx 


\ , (58) 




Jo 





where the disorder average of the exponential involves 
functional integration over V (x). As an example, we 
calculate A sc for Gaussian disorder in the limit |e| <C 
Vo and ^/VqR c 3> 1, which justifies the stationary point 
approximation. The latter yields 



A sc = - (2R C V 2 ) 



1/3 



(59) 



According to Eqs. (|59")l and (|D3|) . under the conditions 
|e| <C Vo and \/VqR c ^> 1, LE ratio p — A/ci scales 
like (R C V 2 ) 1/3 /V 1/2 ~ ^ 1/3 



(cii? c ) J 



For comparison, the 
Gaussian part of p, i.e. = \ ( 1 + fj) [cf. Eq.([T6"|)]. 

scales in this regime as C2/C1 ~ c\R c . The obtained 
scaling is closely related to the statistics of the disorder 
fluctuations and, thus, is specific to the considered Gaus- 
sian models (e.g., for strong binary disorder one expects 
p ~ const). 

Let us note that A sc in Eq. (fB"9"]) matches very closely 
the strong disorder limit of the small- R c approximation 
for A (Sec. IIII C| last paragraph). The latter reads 1 < 

A/Ai (.g) < 1.16, where g = 2V^R C and Ai (3) = \ (f ) 1/3 
is the generalized LE for (5-correlated disorder in strong 
disorder limit (cf. Eq. (|A4|) V The semiclassical result (|59|) 
gives A sc /Ai (g) = 1.19, which is remarkably close to the 
small- R c approximation. Note that the considered limit, 
VVo-Rc 3> 1, implies that AR C ^> 1, which is formally 
beyond validity of the small- R c approximation. 



V. TRANSMISSION STATISTICS 

In this section we discuss properties of the transmis- 
sion coefficient distribution in terms of the dimensionless 
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FIG. 3: (color online) LE ratio p = A/ci as a function of the 
dimensionless disorder intensity /3 [Eq. (|49[) ]. LE ci is ob- 
tained by numerical simulations. GLE A is calculated using 
both the analytical RWNA method, panel (a), as well as nu- 
merical simulations, panel (b). Legends in (a) apply to both 
panels and specify type of correlation and value of kR c . The 
analytical and the numerical results complement each other 
in different regimes of disorder. The numerics is noisy or 
unavailable at sufficiently large disorder intensities, /3 > 0.1, 
where analytical results demonstrate how p increases with 
kR c . On the other hand, numerical results are more accurate 
for smaller disorder intensities, j3 < 0.1, where deviation from 
the limiting weak disorder value p — 1 begins to develop, as 
shown in the inset of panel (b). 



ratios c^/ci, C3/C1 and p = A/ci, Eqs. (fT3| and (|16jl . 
The cumulant coefficients are simulated numerically for 
three types of correlations listed in Table fl] and for the 
white noise model as well. For the generalized LE A 
we use both analytical and numerical results. Details 
of the numerical simulations are given in Appendix [C] 
All calculations are performed at the same fixed value of 
energy, e = k 2 > 0, and for different values of the corre- 
lation radius R c and the disorder variance V Q 2 . Disorder 
strength is conveniently characterized by the dimension- 
less disorder intensity fi = V^kRc/Ae 2 , introduced in Eq. 
(|49|) . For each considered value of R c , quantities A and 
c„ are studied as a function of Vq in the range from weak 
(/3 <C 1) to strong fj3 ~ 1) disorder. 
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FIG. 4: (color online) Cumulant ratios cijc\ (upper panel) 
and cs/ci (lower panel) for the 5-correlated and the expo- 
nentially corelated disorder plotted as a function of j3 (other 
models of correlation given in Table|T]exhibit similar behavior, 
cf. Fig. [6]). Legends in (a) apply to both panels and specify 
the types of correlation and the values of kR c . 



In Figs. [3] and H] we plot p and cumulant ratios c-ijci 
and cs/ci respectively as a function of f3. Different data 
series correspond to different models of correlation and 
different fixed values of kR c , as indicated in the legends. 
All three quantities, p, c^jcx and c^/ci, exhibit quali- 
tatively similar behavior in weak and moderate disorder 
(/3 < 1). As expected from the results on the 5-correlated 
disorder [1],[5|], universal relations p — 1, C2/C1 = 1 and 
cz/ci — (cf. Eqs.dU and (IToTl ) are violated beyond 
weak disorder. Our calculations show that deviation from 
these weak disorder values is strongly enhanced in the 
presence of correlations and increases with the correla- 
tion radius R c . Namely, while this deviation is small 
compared to unity for (5-correlated disorder, it becomes 
of the order of unity, or even larger, in the presence of 
correlations. Below we discuss the obtained results and 
the corresponding parametric dependences. 



A. Scaling of the cumulant ratios c„/ci 

In sufficiently strong disorder, such that ciR c » 1, 
cumulant ratios c n jc\ are described by the semiclassi- 
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FIG. 5: (color online) Semiclassical scaling relations ()57[) are 
verified for the cumulant coeffcients C2 (upper panel) and C3 
(lower panel) by plotting simulated values of c n /c1Rin — 1) 
as a function of the parameter e/Vb- Legends in (b) apply to 
both panels. The dashed lines represent analytical result due 
to Eq. |g7]| with B 2 = 0.8 and B :i = 1.2 (see text after Eq. 
(|57p ). To check the semiclassical approximation for LE ci, 
Eq. (|51l) . inset in panel (a) shows that ratio of the simulated 
value of Ci to its semiclassical approximation, denoted as c\ c , 
approaches unity for e/Vo < 1. As expected, semiclassical 
approximation improves for larger values of R c and smaller 
values of e/Vo- 



cal relations (|56|) . To verify these relations, in Fig. [5] 
we plot (c„/ci) / (ciR c ) n for n = 2,3 as a function 
of e/Vo aid compare it to the semiclassical prediction 

S n -iB n (^^j fn (po) ' accor ding to Eq. (f56|). For Gaus- 
sian disorder, functions /„ {^y^j ar e given by Eq. (|D4I) . 
and we use values B 2 = 0.8 and B 3 = 1.2, as explained 
in text after Eq. (|57|) . Fig. [5] shows that the obtained 
numerical values converge to the theoretical curve for 
e/Vb < 1. As expected, the semiclassical approximation 
improves for smaller e/Vo and larger R c . Let us note, 
for clarity, that in our presentation we increase the dis- 
order amplitude Vq keeping fixed e > and fixed R c . In 
this parametrization, cumulant ratios c n /c\ grow mono- 
tonically with j3, roug hly as (cii? c ) n_i (cf. Fig. i}. If, 
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FIG. 6: (color online) Demonstration of the empirical scal- 
ing (|60l) for the cumulant coeffcients ci (upper panel) and cz 
(lower panel). Legends in (b) apply to both panels and spec- 
ify type of correlation and value of kR c . We use a — for 
the exponential correlation, and a = 0.4 - for the Gaussian 
and the speckle correlations. For better visibility, data for 
the exponential correlation are offset from zero by a constant 
(in both panels). Inset: same as in the main panel (a), but 
plotted in a logarithmic scale for j/-axis and without the offset. 



is not universal as well, and we find that a w 0, for 
the exponential, and a s» 0.4 for the Gaussian and the 
speckle correlation functions (Table [TJ) . On the con- 
trary, in the considered Gaussian models, we obtain 
V2 ~ 3/4 and ^3 ~ 9/8, irrespectively of the disorder 
correlation. In Fig. [5] we plot (02/01 — 1) / (kR c ) 3 ^ 4 

and (cs/ci) / (kR c ) 9 ^ 8 against (kR c )~ a lg c\R c , which 
demonstrates that scaling (ISUl) holds for a rather broad 
range of the considered values of kR c (see legends) 
as long as c\R c < 1. A logarithmic scale plot of 
(ca/ci — 1) / (fci? c ) 3 ^ 4 , given in the inset of the upper 
panel, shows that this scaling holds also at small values 
of c\R c (the same is true also for c^/ci). 

Let us stress that relations (|60)) become meaningless in 
the white noise limit kR c — > 0, and are applicable only 
for kR c > 1, when deviation from the weak disorder rela- 
tions ([4]) is dominated by the effects of the correlations. 
According to Eq. ([5D|) . the latter is controlled by the 
ratio of the correlation to the localization lengths, C\R C . 

As seen in Fig. ©, functions fi n are rather similar, 
though not identical, for different models of correlation. 
The major distinction in a specific form of relation (|60[) 
for different (Gaussian) models appears in the value of 
the exponent a. It is interesting to note that both \i n 
and a practically coincide for the models with Gaussian 
and speckle correlation functions. This can be related to 
the qualitative similarity of their power spectra, given by 
the Gaussian and the "tent" function respectively, which 
have either effective or exact cutoff of the order of R^ 1 
(the "tent" function, (1 — \x\) 9 (1 — \x\), was introduced 
in Sec. IIIIEp . On the contrary power spectrum of the 
exponentially correlated disorder is given by the slowly 
decaying Lorentzian function. Note that disorder power 
spectrum appears, e.g., in the Born approximation (|B3p 
for A (and LE ci), which shows close relation between 
the properties of localization and those of the disorder 
power spectrum. 



B. Extreme fluctuation of T 



however, relative disorder strength is increased by de- 
creasing e at fixed Vq and R c , then ratios c n jc\ would 
eventually vanish for sufficiently large and negative e, as 
can be seen from Eq. (|5o) [because Vq/ \e\ becomes small 
in this limit]. 

Semiclassical scaling (|56[) is not applicable in weak and 
moderate disorder. In this regime, and for kR c > 1, we 
find empirically that deviation from the weak disorder 
relations @ for n — 2,3 is described by the following 
approximate scaling 

t-^=<^Mw)' (60) 

where S n 2 is the Kronecker delta. Here /i„ (x) are some 
dimensionless monotonically increasing functions, whose 
form depend on the type of correlation. Exponent a 



Beyond weak disorder, LE ratio p increases with the 
correlation parameter kR c and, as a function of the dis- 
order intensity, exhibits a peak [42| at moderate disorder 
(Fig. [3]). Thus, as expressed in terms of the ratio between 
the mean and the typical values of the inverse transmis- 
sion coefficient T _1 , Eq. (|17p. also the extreme relative 
fluctuation of T~ l is peaked at moderate disorder and is 
enhanced by the disorder correlations. This conclusion 
is consistent with our observation that statistical conver- 
gence of the Monte-Carlo simulation of the generalized 
LE A is most slow at moderate disorder and for larger 
values of R c , as indicated by the error bars in Figs. Q] 
and[U Note that relative fluctuation (fTF]) depends expo- 
nentially on pi, where I = c\L. Therefore, for kR c > 1 
and I > 1, correlations lead to the exponentially large 
enhancement of the extreme fluctuation of T -1 as com- 
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FIG. 7: (color online) The non-Gaussian part of p, ApNG = 
p — pa, where pa = (1 + c%/ci) /2, as a function of /3. The 
shown data correspond to the 5-correlated and exponentially 
correlated disorder with different values of kR c , as indicated 
in the legends. A is calculated analytically using the RWNA 
method. As noted in Sec. IIII El approximation error of this 
method becomes significant for kR c > 1. In the presented 
figure, this does results in false values for /? < 0.1 (simulated 
values are much closer to zero). However, the approxiamtion 
error is not critical in the region of main interest, 0.1 < /3 < 1, 
where it tends to underestimate, rather than overestimate, the 
actual values of ApNG- Thus, the presented effect is not an 
artifact of the approximation. 



pared to the case of the (^-correlated disorder. 

It is instructive to compare p to its Gaussian part pq = 
| (l + 2a. J. The latter is obtained by discarding the non- 
Gaussian terms, n > 3, in Eq. (|16l) . which is equivalent to 
approximating the asymptotic distribution of In T by the 
Gaussian one with the same mean and variance, equal 
to 2c\L and 4c2^ respectively. Correspondingly, con- 
tribution of the non-Gaussian corrections is represented 
by the difference Ap NG = p - p G = \ £^=3 frff- In 
Fig. [7] we plot AppjQ for the ^-correlated and the expo- 
nentially correlated disorder (other models of correlation 
exhibit similar behavior). In agreement with relations 
(j4|), ApNG vanishes in the weak disorder limit. Beyond 
this regime, Ap^c is positively peaked at moderate dis- 
order, and the magnitude of the peak increases with kR c . 
Positive values of Apj^c mean that the low-T tail of the 
In T-distribution is heavier than in the Gaussian approxi- 
mation in a sense that (T _1 ) and, thus, the relative fluc- 
tuation (|17[) . are larger than in the corresponding Gaus- 
sian distribution. The result in Fig. [7] shows that this 
"super-Gaussian" effect is enhanced by the disorder cor- 
relations and is most prominent at moderate disorder. 
Let us stress, however, that this statement applies only 
to the integral quantity (|17p and does not imply any spe- 
cific form (e.g., super- or sub-Gaussian) for the asymp- 
totic decay law of the low-T tail in the In T-distribution. 
The latter was studied , e.g., in Refs. Q, (43[, 44 1 using 



different methods. 



VI. CONCLUSIONS 

We have considered statistical properties of the trans- 
mission coefficient T of a one-dimensional disordered sys- 
tem described by the Schrodinger Eq. ([T]) with a Gaus- 
sian correlated disorder. The main focus of our study 
was an effect of the correlations on the transmission dis- 
tribution, which was characterized in terms of the di- 
mensionless ratio between the generalized and the usual 
Lyapunov exponents, p = A/ci, as well as ratios cijc\ 
and cs/ci of the asymptotic cumulants of InT -1 . Both 
analytical and numerical methods were employed for cal- 
culation of A and c n . 

First, a small- R c approximation was developed for the 
generalized LE A, which is not limited to weak disorder 
and also is able to account for the non-trivial effects of 
the correlations. To this end, we obtained an infinite 
hierarchy of integro-differential equations on the second 
moments of the wave function and its functional deriva- 
tives with respect to the disorder. We have shown that 
this hierarchical chain can be truncated starting from the 
second level to obtain a non-trivial approximation, which 
accounts for the correlations of the disorder. Note that 
termination of this hierarchy on the first level accounts 
only for the intensity of the noise, while all information 
on the correlation properties is lost. This non-trivial clo- 
sure allowed us to obtain an analytical solution for the 
generalized LE in an implicit form as a largest real root 
of a non-linear algebraic equation. The obtained approx- 
imation is valid formally for kR c <C 1 and AR C <C 1, 
where k = y/\e~\ is the wave number. It turned out that a 
standard white noise approximation (36j | does not allow to 
treat the generalized LE in the case of sub-exponentially 
decaying correlations. To overcome this obstacle, we have 
proposed a simple self-consistent regularization, which 
extends the applicability of the approximation to any 
correlation function satisfying | J Q C 2 (x)dx\ < oo. The 
derived small- R c approximation was compared to numer- 
ical simulations and a good agreement was found. 

The asymptotic cumulants of In T were calculated an- 
alytically within the semiclassical approximation appli- 
cable in a special case of strong disorder, e < Vo and 
RcVVo — e 1. For an arbitrary strength of disorder, 
cumulant coefficients 01,2,3 were obtained by numerical 
simulations. 

Using these analytical and numerical methods, the gen- 
eralized LE and the first three cumulants of InT were 
calculated for several models of correlations and for dif- 
ferent values of R c , which enabled us to investigate effects 
of correlations on the form of the transmission distribu- 
tion. In order to study transition between the regimes of 
weak and strong disorder, we have considered only posi- 
tive energies e, where it was convenient to introduce the 
dimensionless disorder intensity /3, defined in Eq. (|49l) . 

In sufficiently weak disorder, we obtained c%jc\ = 1, 
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ca/cx = and p = 1 in all the considered cases. 
Thus, as expected from previous work (e.g., Refs. 
[^,[^,[23, [2^{), in weak disorder, correlations do not 
destroy SPS and do not modify the universal relations 
((4]). In the white noise limit, kR c -C 1, this regime is re- 
alized for ci/fc ~ j3 <C 1. For kR c > 1, the relevant con- 
trol parameter is the ratio of the correlation to the local- 
ization lengths, c\R c . Quantity (cii? c )~ 1 is interpreted 
naturally as a measure of randomization of the disorder 
potential on the scale of the localization length. Corre- 
spondingly, the weak disorder universality, expressed by 
relations (UJ), takes place when C\R C <C 1 (for kR c > 1, 
this means that disorder is weak also in the conventional 
sense ci/fe -C 1). 

Relations (HJ) are not valid beyond weak disorder. In 
the white noise model, the corresponding deviation from 
these relations is weak for any strength of the disor- 
der (for e > 0), in a sense that values of p and Cijc\ 
stay very close to unity, while C3/C1 remains small com- 
pared to unity. On the contrary, in correlated disorder, 
these quantities depend strongly on the strength of the 
disorder through the control parameter c\R c . Regard- 
ing this parametric dependence, we discuss two regimes: 
c\R c < 1 and ciR c ^> 1. 

The first regime, c\R c % 1, corresponds to transition 
from weak to moderate disorder. In this regime and for 
kR c ^ 1, cumulant ratios cijc\ and c^/ci are described 
by the approximate scaling relations (|60[) with param- 
eters kR c and C\R C . This relations demonstrate that, 
starting from weak disorder (c\R c <C 1), ratios c n /ci in- 
crease gradually with parameter c\R c and, for c\R c ~ 1, 
arrive at values which can be much larger than unity 
(depending on kR c ). In such a case, the bulk of the lnT- 
distribution becomes much broader than in weak or in 
the white noise disorder. 

While the small-n cumulant coefficients c n describe 
bulk of the In T-distribution, LE ratio p is a measure 
of the extreme relative fluctuation of T _1 , expressed by 
the ratio between the mean and the typical values of T _1 , 
Eq. (|17p . Like the cumulant ratios, p increases from weak 
to moderate disorder and its deviation from the weak 
disorder value p = 1 is strongly enhanced by the dis- 
order correlations (namely, exceeds unity for kR c > 1). 
As a function of the disorder strength, p and, thus, the 
extreme relative fluctuation (|17[) , are peaked at moder- 
ate disorder (near c\R c ~ 0.1). This peak of p is as- 
sociated with the non-Gaussian corrections to the low- 
T tail of the In T-distribution, whose contribution to p 
is positively peaked, i.e. "super-Gaussian", at moder- 
ate disorder. The latter has the following simple inter- 
pretation in terms of the disorder statistics. In moder- 
ate disorder, when energy e is of the order of Vo, wave 
propagation becomes affected by the under-barrier tun- 
neling through the rare but large peaks of the random 
potential. For ^/V^R C ~ kR c > 1, already a single large 
barrier becomes a strong scatterer. Therefore, fluctua- 
tions in height and in occurrence of these rare peaks lead 
to extreme deviations of T -1 from its typical value. In 



stronger disorder, the typical value grows significantly, 
and relative contribution of the large rare barriers be- 
comes less pronounced. This simple explanation can be 
confirmed by calculating p in the framework of the semi- 
classical approximation (|50|) [not presented here] . 

The second regime, c\R c 3> 1, is realized when e < Vo 
and R c y/Vo — e 3> 1. In these conditions, interference 
effects are suppressed and localization is dominated by 
under-barrier tunneling. Then, according to Eq. (|50l) . 
the transmission distribution is directly related to the 
statistics of the "excursions" of the random potential 
above the level V (x) — e, and the cumulants of lnT -1 
satisfy the " semiclassical" scaling relations (l56|) , which 
were verified in numerical simulations. According to Eq. 
([56| . the effect of correlations is expressed by the sim- 
ple relation c n jc\ ~ f n ( e /Vo) (ciR c ) n 1 , where the co- 
efficient /„ (e/Vb) = ((k™)) / depends only on the 
one-point distribution of disorder. Semiclassical approx- 
imation can also be used to calculate the generalized LE 
A. In particular, for |e| -C Vq and y/Vo~R c ^> 1, we ob- 



tain scaling p ~ {c\R L 
statistcs of disorder. 



a/3 



, which is specific to Gaussian 
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Appendix A: Exact solution for A in 5-correlated 
disorder 



Equation (1251) decouples and can be solved exactly for 
the (5-correlated disorder (see e.g. Ref. Q). Namely, sub- 
stituting O2 (t — t) = g8 (t — t) into ([23]) and using the 



initial condition 



SY(t 



sv(t-) 

d t (Y(t))=M(e,g) (Y (i)) , 



V (Y (i)) (27]), one obtains 
(Al) 



where 



M(e,g)=C + ^V 2 = V2 



1 

-e 1 

g/V2 -e 



(A2) 



and matrices C and T> are defined in (|2Tj) . Solution of 
Eq. Ig]} is 



(Y (*)) = e^'^Yo, 



(A3) 



where Yo = Y (0) is an initial condition. The generalized 
LE A, Eq. (fl4l) . is given by the largest real eigenvalue of 
the matrix (e, g). It is found from the cubic equation 



A(4A 2 



9 



(A4) 
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which has the following roots 



Ai 



3p {g, e) ! 



A 



2.3 



-Ai ± V 3A i + 6 



where p (5, e) = (g + \J ' g 2 + (§e) 



1/3 



(A5) 

The general- 



ized LE A is given by Ai, which is real and positive. 
Other two eigenvalues are complex for e > — 3A 2 (i.e. 
e > — fg 2 / 3 ), and real otherwise. Matrix M. (e, g) is not 
normal (M.Mr 7^ JvftM), and its left and right eigen- 
vectors, corresponding to the eigenvalues Aj in Eq. 
are 



5|- 



,R 



1 



24A? + 2e 



1, 2y/2ki, e + 8A? 



e + 8Af, 2^Ai, 1 



(A6) 



These eigenvectors satisfy the normalization (w^itj 1 
5u- 



Appendix B: Born approximation for A 

The systematic weak disorder expansion (i.e. in pow- 
ers of the disorder amplitude Vo) for the generalized LE 
A was considered in Ref . [22[ . Here we only note that the 
Born approximation for A can be obtained by substitu- 
tion of the zero order solution for ^ |^^y ^ into (|23l) . The 

zero order solution for / \ is obtained by neglecting 
the noise term in equation (|26p : 

9 /SY(t)\ =c /SY(t)\ / SY(t+) \_ 



dt\sV(T)/ \sV{t)/' \ 6V(t) 
which has the solution 



= ©(Y(r)), 



(Bl) 



Inserting this solution into Eq. (|23|) yields an equation 
closed with respect to (Y (t)): 

8 t (Y) = C (Y) + [ C 2 (t - t) Ve^-^V (Y (r)) dr, 
Jo 

(B2) 

where pair correlation function C2 (t — r) is defined in 
Eq. (J6)). Then, lowest order calculation of the asymp- 
totic growth rate of the solution (Y (t)) yields the Born 
approximation 



A 9 



k 8e 3 / 2 



J dsT{s) cos {2kR cS ) + o(-^-j , (B3) 



where k = *Je and g = 2V 2 R C . This expression coincides 
with the Born approximation for LE c\ Q • 



Appendix C: Method of numerical simulations 

The numerical simulations were performed using the 
tight-binding (TB) model 



Eip n = e n ip n - ip n+1 - tpn-i, 



(CI) 



with energy near the band edge, < E + 2 <C 1, and 
the diagonal disorder \e n \ < E + 2. In this regime, Eq. 
(jCljl is a good approximation to the continuous model 
([5]). We relate the continuous model with energy e = k 2 
and potential V (x) to the TB counterpart by setting 



E = — 2cos/c, 



sin k 



V(x)\ 3 



(C2) 



The correlated Gaussian disorder e„ was generated by 
filtering sequences of independent Gaussian random vari- 
ables. The real space convolution with a proper ker- 
nel was used to obtain the "short-range" exponential 
(r (x) = e~ x ) and Gaussian (r [x) — e~ 7TX / 4 ) correla- 
tions, while the Fourier space filter was applied to obtain 
the power law correlation T (x) = sine 2 | . 

In all simulations presented in this paper, we have fixed 
the energy at E = —1.95 and varied values of the disorder 
amplitude Vq and the correlation radius R c . Note that for 
E = —1.95, the corresponding wavelength of the solution 
for a pure system (e„ = 0) is equal to about 28 sites. 
Therefore, approximation to a continuum is good as long 
as Vo < E. 

A standard transfer matrix formalism (see e.g. Ref. 
[45j |) was used, in which the TB equation (|C1|) is rewrit- 
ten in the matrix form 




-T ^ n I T - 

— In I 1 I j In — 

\ Wn-l J 



(e n -E) -1 
1 



(C3) 

where 7~ n is a single-site transfer matrix. Then, solution 
of the initial value problem is given by 




N,l 




Tn,i=T N -T 2 Ti, (C4) 



where 7i,jv is the total transfer matrix for the system of 
length N. 

In analogy with Eq. (|19p . one can introduce vector 
Y„ = (u 2 , \/2u n v n , w 2 ) T , where u» = ip n and v n = 
(tp n — ip n -i) /k ~ k~ 1 d x i(). Similarly to (|C4|) . the so- 
lution for Yjy can be write as 



Yat+i — /Cjv.iYi, 



(C5) 



where Yi is an initial condition, and /Cjv.i is the transfer 
matrix for Y^v , which is readily expressed in terms of the 
elements of 7jv,i- The largest eigenvalue of ICn, 1 is equal 
to the square modulus of the largest eigenvalue of Tn,i- 
Finally, ensemble average over the disorder realizations 
yields 



(Y 



N+l/ 



(ICn^Yl 
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Cumulant coefficients c n , denned in Eq. (fl"2|) . are given 
by the asymptotic linear growth rate of the cumulants 
2~ n ((In" Kjv)) with the system length N, where kn de- 
notes the largest eigenvalue of the matrix 1C\^. The 
asymptotic slope was calculated by the linear fit, which 
have to exclude the region of the initial transient of the 
order of a few localization lengths. The ensemble average 
was performed over ~ 10 6 -j- 10 7 realizations of disorder. 

The generalized Lyapunov exponent A, Eq. (|14p. was 
calculated as a linear slope of j In (kn) ■ Alternatively, A 
could be found as a slope of the logarithm of the largest 
eigenvalue of (/Cjv.i), which gives practically the same 
result. About 10 s realization were generated to calculate 
each value of (kn)- 

Monte- Carlo simulation of the generalized LE can be 
a quite challenging task, as is briefly explained in the 
following. In numerical simulation of the generalized LE 
one have to deal with two restrictions on the system size, 
both from below and from above. The lower bound is de- 
termined by the width of the transient to the asymptotic 
behavior 



tually coincide in sufficiently strong disorder, since c^ 1 
becomes small. In such a case, numerical calculation of A 
becomes impossible, unless the number of the realizations 
is increased dramatically (for the exactly log-normal kn, 
it can be shown that the improvement is logarithmically 
slow). In moderate disorder, the small range between the 
lower and the upper bounds results in large uncertainty 
in the calculated A, as indicated by the error bars in Figs. 
Q]and[2j 



Appendix D: Auxiliary formulas for semiclassical 
approximation 

Cumulants ((«")) are easily written in terms of ((«;") ) e 
and ((9™ {V))) using that 9™ (V) = 9 e (V) and (« n ) = 
(k™} ( (9 e (V)), where (• • • ) stands for the average with 
the distribution P e (V) = 9 e (V) P (V) / / e °° dVP (V), 
introduced in Sec. IIVI For example, the second cumu- 
lant, ((k 2 )) = (k 2 ) — (k) 2 , becomes 



In (kn) = NA + const. 



(C6) 



Width of this transient is at least of the order of R c . This 
follows from the form of the differential equation (|23[) for 
(Y (t)), which suggests that growth rate of (Y (£)) can not 
stabilize unless t 3> R c - This is because the correlation 
function in the integral on the right hand side of Eq. 
(l23l) decays on the scale of R c . In the white noise limit, 
R c — > 0, the transient region is absent, as follows from the 
exact analytical result (Appendix [X]) and was observed 
numerically. Therefore, we assume that the width of the 
transient is of the order of R c , and other scales, such 
as the localization length, are less important (unlike the 
case of Cn, n > 2). 

The upper bound on the system size is determined em- 
pirically from the numerical data as a value of N, beyond 
which the linear dependence in Eq. (|C6j) is violated. 
This computational artifact originates from the insuffi- 
cient statistics in averaging of the broadly distributed 
quantity /cjv, which is equivalent to T _1 . According to 
Eqs. ((2} and (|4]), distribution of njy is log- normal in 
weak disorder, while some corrections to the limiting log- 
normal form appear in stronger disorder. As follows from 
Eqs. (fl3|) and (|17|) . this distribution becomes increas- 
ingly broad and heavy-tailed with the increase of the di- 
mensionless system length I = c\N . For large I, long tails 
of the distribution, which dominate the theoretical mean 
of /-cat, are typically under-sampled in simulations with 
a finite number of realization. As a result, the obtained 
values of In (kn) become typically underestimated (for- 
mally, expectation value of In (kjv) becomes smaller that 
ki^K/v)^ where (hn)^ is the theoretical mean). This 
effect increases with the system size, the relevant length 
scale being the localization length c^ 1 . Therefore, the 
upper bound on the system length is of the order of a 
few c^f 1 . Thus, the upper and the lower bounds even- 



«k 2 )) = < K 2 ) £ (9 £ (V))-(n) 2 e (e e {V)Y 



(Dl) 



which, using ((9 2 (y))) = (e e (V)) - (9 e {V)) 2 and 
(( k2 )) c = ( k2 ) t~ yi^ds after some rearrangement 

((- 2 )) = ((- 2 ))A^(V)) + (n) 2 m(V)))- (D2) 

Similar expansion for the third cumulant, ((k 3 )) = 
(k 3 ) — 3 (k 2 ) (k) + 2 (k) 3 , gives expression (|57|) for f 3 — 
«- 3 ))/(-> 3 . 

For Gaussian distribution P (V) = 

(2ti% 2 )~ 1/2 e~ v2 / 2V ° , definition (J51J) yields the fol- 
lowing expression for LE: 



d = Vj /2 2-^e- r2 / 2 



u 



3 1 r 2 
4' 2' ~2 



(D3) 



while functions /„ = / f° r n — 2,3 are 

given by (r > 0) 



h(r) = -l 



2e r2 / 2 - V2nre r2 erfc (r/V2) 



h(r) 



2e— 



24- 



2irre r erfc -j= 



+ 



32rae 4 



(D4) 



where U (a, b, z) is a confluent hypergeometric function 
of the second kind (or Tricomi function) and K v (z) is a 
modified Bcssel function of the second kind I4a| . 
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